PRL  101,  057003  (2008) 


PHYSICAL  REVIEW  LETTERS 


week  ending 
1  AUGUST  2008 


Unconventional  Superconductivity  with  a  Sign  Reversal 
in  the  Order  Parameter  of  LaFeAsO^F* 

1. 1.  Mazin,1  D.  J.  Singh,2  M.D.  Johannes,1  and  M.H.  Du2 

1Code  6393,  Naval  Research  Laboratory,  Washington,  D.C.  20375,  USA 
2 Materials  Science  and  Technology  Division,  Oak  Ridge  National  Laboratory,  Oak  Ridge,  Tennessee  37831-6114,  USA 

(Received  18  March  2008;  published  29  July  2008) 

We  argue  that  the  newly  discovered  superconductivity  in  a  nearly  magnetic,  Fe-based  layered 
compound  is  unconventional  and  mediated  by  antiferromagnetic  spin  fluctuations,  though  different 
from  the  usual  superexchange  and  specific  to  this  compound.  This  resulting  state  is  an  example  of 
extended  s-wave  pairing  with  a  sign  reversal  of  the  order  parameter  between  different  Fermi  surface 
sheets.  The  main  role  of  doping  in  this  scenario  is  to  lower  the  density  of  states  and  suppress  the  pair¬ 
breaking  ferromagnetic  fluctuations. 
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The  discovery  of  superconductivity  with  Tc  S  26  K  [  1  ] 
in  a  compound  that  contains  doped  Fe2+  square  lattice 
sheets  raises  immediate  questions  about  the  nature  of  the 
superconducting  state  and  the  pairing  mechanism.  A  num¬ 
ber  of  highly  unusual  properties  suggest,  even  at  this  early 
stage,  that  the  superconductivity  is  not  conventional.  We 
argue  that  not  only  is  it  unconventional  but  that  doped 
LaFeAsO  represents  the  first  example  of  multigap  super¬ 
conductivity  with  a  discontinuous  sign  change  of  the  order 
parameter  (OP)  phase  between  bands,  a  state  discussed 
previously  (e.g.,  Refs.  [2,3])  but  not  yet  observed  in  nature. 

We  suggest  that  superconductivity  here  is  mediated  by 
spin  fluctuations  (SFs),  as  many  believe  is  the  case  in 
cuprates,  heavy  fermion  materials,  or  ruthenates.  SFs  can 
induce  only  a  triplet  superconducting  OP  or  a  singlet  one 
that  changes  sign  over  the  Fermi  surface  (FS).  The  latter 
condition  often,  but  not  always,  dictates  strong  angular 
anisotropy  of  the  OP  (cf.  d  wave).  In  our  scenario,  it  is 
satisfied  despite  full  angular  isotropy,  since  the  sign 
changes  between  the  two  sets  of  FSs.  Our  model  is  also 
principally  different  from  conventional  .y-wave  supercon¬ 
ductivity  discovered  in  MgB2:  There  the  pairing  interac¬ 
tion  is  attractive;  in  our  case  it  is  repulsive  (but  pairing 
thanks  to  the  sign  reversal).  Finally,  similar  to  d-  or  p-wave 
superconductivity,  our  OP  has  a  nearest-neighbor  structure 
in  real  space,  thus  reducing  the  Coulomb  repulsion  within 
the  pair. 

We  begin  by  arguing  against  conventional  superconduc¬ 
tivity.  The  pure  compound  LaFeAsO  is  on  the  verge  of  a 
magnetic  instability:  It  has  a  very  high  magnetic  suscepti¬ 
bility  [1]  and  is  strongly  renormalized  compared  to  density 
functional  (DFT)  calculations  [4].  This  renormalization  is 
higher  than  in  any  known  conventional  superconductor, 
including  MgCNi3,  where  superconductivity  is  believed 
to  be  strongly  depressed  by  spin  fluctuations.  The  suscep¬ 
tibility  in  the  pure  compound  is  large  and  upon  doping  with 
F  grows  even  larger  and  becomes  Curie -Weiss-like.  This 
suggests  nearness  to  a  critical  point  in  the  pure  compound 
and  nontrivial  competition  between  different  SFs.  Very 
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strong  electron-phonon  interactions  would  be  required  to 
overcome  the  destructive  effects  of  such  SFs.  We  have 
calculated  ab  initio  the  electron-phonon  spectral  func¬ 
tion  a2F(co)  and  coupling  A  for  the  stoichiometric  com¬ 
pound  [5].  Some  moderate  coupling  exists,  mostly  to  As 
modes,  but  the  total  A  appears  to  be  —0.2,  with  <ulog  — 
250  K  (rui0g  is  the  logarithmically  averaged  boson  fre¬ 
quency  of  Eliashberg  theory),  which  can  in  no  way  explain 
Tc  >  26  K. 

The  calculated  DFT  Fermi  surfaces  [4]  for  undoped 
LaFeAsO  consist  of  two  small  electron  cylinders  around 
the  tetragonal  M  point  and  two  hole  cylinders,  plus  a  heavy 
3D  hole  pocket  around  I  .  To  study  doping  effects,  we 
performed  full-potential  calculations  using  the  WIEN  pack¬ 
age  in  the  virtual  crystal  approximation  and  the  Perdew- 
Burke-Ernzerhof  generalized  gradient  approximation 
(GGA)  functional  [6].  The  lattice  parameters  we  took 
from  experiment  [1]  and  optimized  the  internal  positions 
[7].  In  Fig.  1,  we  show  the  bands  near  the  Fermi  level  for 
x  =  0. 1  and  the  corresponding  Fermi  surface.  As  expected, 
the  3D  pocket  fills  with  electron  doping  (at  x  =  0.04—0.05) 
and  the  fermiology  radically  simplifies,  leaving  a  highly 


FIG.  1  (color  online),  (a)  Calculated  band  structure  at  x  =  0. 1 
near  the  Fermi  level,  (b)  Calculated  Fermi  surface  at  10% 
doping.  Note  that  the  only  3D  parts  are  the  far  ends  of  the 
electron  cylinders  around  M.  The  fully  three-dimensional  sur¬ 
face  present  in  the  undoped  compound  is  suppressed  beneath  EF 
by  increased  electron  count. 
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2D  electronic  structure  with  two  heavy  hole  cylinders  and 
two  lighter  (and  larger)  electron  cylinders. 

This  fermiology  imposes  strong  constraints  on  the 
superconductivity.  In  particular,  with  the  exception  of  pho¬ 
nons,  it  is  hard  to  identify  pairing  interactions  with  a  strong 
kz  dependence.  Thus,  states  with  strong  variations  of  the 
OP  along  kz  are  unlikely.  An  angular  variation  of  the  OP  in 
the  xy  plane  is  possible  but  would  require  an  unrealistically 
strong  q  dependence  of  the  pairing  interaction  on  the  scale 
of  the  small  Fermi  surface  size  and  would  also  be  ex¬ 
tremely  sensitive  to  impurities. 

On  the  other  hand,  the  small  Fermi  surfaces  are  readily 
compatible  with  a  pairing  state  with  weak  variations  of  the 
OP  within  the  sheets  but  a  tt  phase  shift  between  electron 
and  hole  cylinders.  Here  we  show  that  a  SF  pairing  inter¬ 
action  favoring  exactly  such  a  state  is  present  in  this 
material,  and  we  discuss  the  expected  consequent  physical 
properties. 

The  SF  spectrum  is  unusually  rich  for  this  compound 
and  comes  from  three  separate  sources.  First,  the  system 
is  relatively  close  to  a  Stoner  ferromagnetic  instability. 
Second,  there  is  a  nearest-neighbor  antiferromagnetic 
(AFM)  superexchange  [8(d)].  Third,  there  are  nesting- 
related  AFM  spin-density-wave-type  SFs  near  wave  vec¬ 
tors  connecting  the  electron  and  hole  pockets.  The  latter 
appear  to  be  the  strongest  ones.  The  corresponding  inter¬ 
action  connects  the  well-separated  FS  pockets  located 
around  T  and  around  M.  Though  repulsive  in  the  singlet 
channel,  these  would  nevertheless  be  strongly  pairing  pro¬ 
vided  that  the  OPs  on  the  two  sets  of  the  FSs  have  opposite 
signs.  The  main  message  of  our  Letter  is  that  this  “s+” 
superconducting  state  is  both  consistent  with  experimental 
observations  and  most  favored  by  SFs  in  this  system. 

As  opposed  to  the  undoped  material  [4],  we  do  not  find 
any  FM  solution  for  the  doped  compound,  even  with  the 
more  magnetic  GGA  functional.  This  suggests  that  the 
main  function  of  doping  is  to  move  the  system  away 
from  a  ferromagnetic  instability  (see,  however,  Ref.  [8]). 
That  the  system  becomes  less  magnetic  is  due  to  the  fact 
that  upon  doping  the  heavy  3D  hole  pocket  near  T  rap¬ 
idly  fills  and  the  total  density  of  states  (DOS)  drops 
(by  a  factor  of  2).  At  a  doping  level  x  =  0.1,  we  obtain 
in  the  GGA  an  unrenormalized  Pauli  susceptibility 
ReYolq  =  0,  co  =  0)  ~  4  X  10-5  emu/mole  [(V(0)  = 
0.64  states/eV/spin/Fe].  Using  1.1  eV  for  the  Stoner  I 
on  Fe,  we  obtain  a  renormalized  ^(0)  of  0. 14  X 
10-3  emu/mole,  a  large  renormalization  but  much  smaller 
than  what  is  needed  to  explain  the  experimental  value  [1]. 
Note  that  in  the  undoped  system  the  calculated  suscepti¬ 
bility  is  larger,  and  the  experimental  one  smaller,  than  in 
the  doped  one.  This  suggests  that  besides  FM  SFs  there  are 
other,  more  important  spin  excitations  in  the  system. 

The  first  candidate  for  these  is  a  superexchange  corre¬ 
sponding  to  the  standard  simple  nearest-neighbor  checker¬ 
board  antiferromagnetism.  Importantly,  there  is  also 
substantial  direct  Fe-Fe  overlap  [4],  which  leads  to  an 


additional  AFM  exchange  of  comparable  strength  and 
with  the  same  checkerboard  geometry. 

Further  discussion  requires  an  understanding  of  the 
fermiology  in  clearer  terms.  The  band  structure  may  at 
first  appear  intractably  complex,  but  it  is  in  fact  relatively 
simple,  involving  only  three  Fe  orbitals  near  the  Fermi 
level.  The  hole  pockets  around  T  originate  from  Fe  dxz 
and  dyz  bands  that  are  degenerate  at  I  and  form  two  nearly 
perfect  concentric  cylinders.  The  electron  surfaces  are 
better  understood  if  we  recall  that  the  underlying  Fe  layer 
forms  a  square  lattice  with  the  period  d  =  o/V 2  (Fig.  2). 
The  corresponding  2D  Brillouin  zone  (BZ)  is  twice  larger 
and  rotated  by  45°.  If  we  could  “unfold”  the  Fermi  surface 
of  Fig.  1 ,  we  would  find  the  same  two  hole  pockets  at  f  and 
one  electron  pocket  around  the  point  X  in  the  large  BZ. 
The  latter  is  formed  by  the  dyz  band  (or  dxz  near  Y)  that 
starts  180  meV  below  the  Fermi  energy,  disperses  up  along 
X  M ,  and  is  practically  flat  along  X  T .  This  band  is  hybri¬ 
dized  with  the  dxy  band  (or  dx 2  2  in  the  two-Fe  cell), 
which  starts  from  X  at  an  energy  of  —520  meV  and  is 
instead  dispersive  along  X  T .  Upon  hybridization,  these 
two  bands  yield  an  oval  cylindrical  electron  pocket,  elon¬ 
gated  along  f  X  [9]. 

Electronic  transitions  from  the  hole  pockets  to  the  elec¬ 
tron  pocket  should  lead  to  a  broad  (because  the  electron 
pocket  is  oval  rather  than  circular)  peak  in  the  noninteract¬ 
ing  susceptibility  Xo(q>  0J  0),  at  q  =  (77/ a,  0),  while  the 
superexchange  interaction  J{q )  on  the  square  lattice  should 
be  peaked  at  q  —  (7 r/a,  tt/o)  (T  in  the  downfolded 
BZ).  The  renormalized  susceptibility  %{q)  =  ^o(^)/[l  — 
J(q)X o(r?)]  then  has  a  rich  structure  with  maxima  at  f,  X, 
and  M.  For  the  true  unit  cell  with  two  Fe,  both  f  and  M  fold 


FIG.  2  (color  online).  Fermi  surface  formation  upon  backfold- 
ing  of  the  large  BZ  corresponding  to  a  simple  Fe  square  lattice, 
(a)  Real  space:  the  four  small  unit  cells  of  the  Fe-only  sublattice 
(dashed  lines)  with  the  larger  solid  diamond  of  actual  two-Fe 
unit  cell.  Dark  and  light  circles  indicate  superexchange  (checker¬ 
board)  ordering.  The  inset  shows  the  spin-density  wave  corre¬ 
sponding  to  X  point  SFs.  (b)  Reciprocal  space:  the  dashed  black 
square  is  the  unfolded  BZ,  the  solid  blue  diamond  is  the  down- 
folded  zone,  and  the  blue  ellipse  is  the  electron  pocket  from  the 
Y  point  downfolded  onto  the  X  point  (which  is  M  in  the  small 
BZ). 
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down  into  the  T  point,  while  X  folds  down  into  the  M  point. 
The  corresponding  folding  of  the  Fermi  surfaces  makes  the 
electron  pockets  overlap,  forming  two  intersecting  surfaces 
(Fig.  1).  It  is  important  to  appreciate  from  this  gedcinken 
unfolding  that  already  on  the  level  of  the  noninteracting 
susceptibility  there  is  a  tendency  for  antiferromagnetic 
correlations  with  a  wave  vector  different  from  the  super¬ 
exchange  one. 

In  Fig.  3,  we  present  the  calculated  [10]  yo(<7>  w)  = 

at  <°  — '  0  and  ch  =  Vc  (x  is  practically 

independent  of  qz).  The  peak  at  M,  derived  from  interband 
transitions,  is  very  broad,  as  expected,  with  some  structure 
around  the  M  point,  due  to  the  particular  orientation  of  the 
two  oval  pockets  at  M  and  the  size  difference  between  the 
hole  and  electron  cylinders  for  finite  doping.  With  minor 
modification,  this  structure  is  present  also  in  the  undoped 
compound. 

We  have  also  performed  magnetic  calculations  (dis¬ 
cussed  in  a  separate  publication)  in  a  supercell  correspond¬ 
ing  to  the  superexchange  q  =  M  (=  T  in  the  downfolded 
Brillouin  zone)  and  nesting-induced  q  =  X  (=AZ)  spin- 
density  waves.  We  find,  indeed,  that  the  tendency  to  ferro¬ 
magnetic  ordering  is  suppressed  in  the  doped  compound, 
while  the  tendency  to  nesting-based  antiferromagnetism  is 
present,  even  leading  to  an  actual  instability  at  the  mean 
field  level.  It  is  worth  noting  that,  while  this  instability  does 
not  appear  in  actual  superconducting  materials,  it  has  been 
now  found  experimentally  in  the  low  doping  regime  [8]. 

The  strong  AFM  SFs  around  M  favor  our  proposed  in¬ 
state.  Cases  where  SF-induced  interactions  connect  two 
pockets  of  the  Fermi  surfaces,  including  SFs  originating 
from  electronic  transitions  between  the  very  same  pockets 
have  been  considered  in  the  past  [11].  However,  they 
involved  FS  pockets  related  by  symmetry,  which  strongly 
restricts  the  phase  relations  between  them.  In  our  case,  the 
two  sets  of  pockets  are  not  symmetry  related,  and  nothing 
prevents  them  from  assuming  arbitrary  phases. 

For  the  singlet  case,  the  coupling  matrix  between  the 
hole  (h)  and  the  electron  pockets  (e)  is  negative:  Aeh  <  0. 
The  diagonal  components  emerge  from  competition  be- 


FIG.  3  (color  online).  The  imaginary  (a)  and  the  real  (b)  parts 
of  the  noninteracting  susceptibility  XolQ-  w  ~ *  0),  in  arbitrary 
units.  Within  common  approximations  Im^  =  Im[(W(  ^  (  wj] 
is  measurable  by  neutron  scattering  and  Re^  controls  the  pairing 
interaction,  in  the  singlet  channel  proportional  to  1/[1  — 
•f(q)Yo(<7'  0)]  (see  Ref-  [18]  f°r  a  review).  Note  that  the  RPA 
enhancement  will  strengthen  both  peaks. 


tween  the  attractive  phonon-mediated  and  repulsive  SF- 
mediated  interactions  and  are,  presumably,  weak.  If  is 
the  average  of  the  phonon-mediated  interaction  over  the 
wave  vectors  q  <  0  <  2 khF,  and  Xfih  the  same  for  the  SF, 
then  Xhh  =  A^  —  \fh.  Our  calculated  electron-phonon 
coupling  comes  mostly  from  small  wave  vectors,  that  is, 
xfl  +  A?e  ~  0.2  and  A^  ~  0.  Phonons,  therefore,  though 
weak,  promote  the  s+  state.  On  the  other  hand,  Xflh  and  XsJe 
come  from  FM  (small- q)  SFs  and  are  pair-breaking. 
Finally,  the  superexchange  interaction  does  not  affect 
Xhh,  but  the  corresponding  wave  vector  q  =  1  M  connects 
the  e  pockets  near  X  and  Y.  and  such  fluctuations  are  also 
pair-breaking  [8(c)]. 

The  Tc,  as  usual  for  a  two-gap  superconductor  [12],  is 
defined  by  the  maximal  eigenvalue  of  the  A  matrix  Aeff  = 

[A hh  +  Ke  +  \I^XehXhe  +  ( Xee  -  Xhh)2]/2,  and  the  OPs 
A e  h  are  defined  by  the  corresponding  eigenvector.  In  our 
case,  the  signs  of  A/,  and  Ae  will  be  opposite,  and  their 
absolute  values  will  be  similar  since  presumably  (Xee  — 
Xhh)2  Ae/lAAe.  This  “s+”  state  is  analogous  to  the  states 
proposed  previously  for  semimetals  [13]  and  bilayer  cup¬ 
rates  [14]. 

From  the  point  of  view  of  neutron  scattering,  the  struc¬ 
ture  of  the  peak  in  x  near  M  is  important,  but,  for  SF- 
induced  superconductivity,  it  does  not  matter  at  all.  A  well- 
defined  spin  excitation  requires  a  sharp  peak,  but  the  pair¬ 
ing  interaction  is  integrated  over  all  possible  q  vectors 
spanning  the  two  sets  of  Fermi  surfaces  so  that  only  the 
total  weight  of  the  peak  is  important.  Dimensionality  is, 
however,  very  important,  as  3D  coupling  supports  a  long- 
range  magnetic  ordering  (competing  with  superconductiv¬ 
ity).  Indeed,  experiment  shows  an  abrupt  appearance  of 
superconductivity  at  x  ~  0.03,  roughly  where  the  3D  pock¬ 
ets  around  T  disappear. 

One  might  envision  a  triplet  state  similar  to  the  de¬ 
scribed  singlet  one,  fully  gapped,  as  expected  for  unitary 
2D  /j-wave  states,  and  with  different  amplitudes  (or  signs) 
on  the  two  cylinders.  The  similarity,  however,  is  mislead¬ 
ing.  In  the  triplet  channel,  SFs  induce  attraction,  but  given 
the  relatively  large  width  of  the  AF  peak  (Fig.  3),  a  large 
part  of  the  pairing  will  be  lost  as  only  SF  with  a  wave 
vector  exactly  equal  to  (v/a,  n/a)  will  be  fully  pairing, 
and  some  others  will  even  be  pair-breaking.  Therefore,  in 
this  scenario  where  the  antiferromagnetic  SFs  around  the 
M  point  provide  the  primary  pairing  interaction,  we  expect 
the  lowest  energy  superconducting  state  to  be  s+.  The 
structure  of  the  OP  in  real  space  can  be  evaluated  using 
the  lowest  Fourier  component,  compatible  with  the  pro¬ 
posed  s±  state,  namely,  cos kx  +  cosky  (equivalently, 
coskx  cos £,).  This  corresponds  to  pairing  of  electrons  that 
reside  on  the  nearest  neighbors,  just  as  in  the  d- wave  case, 
so  that  the  on  site  Coulomb  repulsion  is  not  particularly 
destructive  for  this  state. 

Finally,  we  discuss  the  experimental  ramifications  of  the 
s±  state.  In  many  aspects  they  are  similar  to  that  of  the  a’± 
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state  proposed  in  Ref.  [14]  and  discussed  in  some  detail  in 
Refs.  [14-16].  The  thermodynamic  and  tunneling  charac¬ 
teristics  are  the  same  as  for  a  conventional  two-gap  super¬ 
conductor  [12].  The  two-gap  character,  however,  may  be 
difficult  to  resolve,  given  the  dominance  of  the  interband 
interactions  which  will  render  the  two  gap  magnitudes 
similar.  Nonmagnetic  small-*/  intraband  (e-e  or  h-h)  scat¬ 
tering  as  well  as  the  interband  spin-flip  pairing  will  not  be 
pair-breaking,  but  paramagnetic  interband  scattering  will, 
resulting  in  a  finite  DOS  below  the  gap  [2],  consistent  with 
specific  heat  measurements  [17].  The  most  interesting 
feature  of  the  s±  state  [16]  is  that  the  coherence  factors 
for  exciting  Bogliubov  quasiparticles  on  FS  sheets  with 
opposite  signs  of  the  OP  are  reversed  compared  to  conven¬ 
tional  coherence  factors.  We  outline  a  few  important  con¬ 
sequences  of  the  relatively  straightforward  application  of 
this  concept  to  experimental  probes.  First,  one  expects  a 
qualitative  difference  between  experiments  that  probe  ver¬ 
tical  transitions  (q  =  0)  and  those  that  probe  transitions 
with  q  close  to  7 t/ a,  i t/ a.  For  instance,  the  spin  suscepti¬ 
bility  at  q  =  0  will  behave  conventionally,  i.e.,  exponen¬ 
tially  decay  below  Tc  without  any  coherence  peak,  while 
the  susceptibility  for  q  ~  77/a,  77/a  will  have  a  coherence 
peak  that  should  be  detectable  by  neutron  scattering  as  an 
enhancement  below  Tc  [16].  AFM  SFs  near  77/a,  77/a 
dominate  in  the  doped  material,  and  so  the  usual  coherence 
peak  in  the  NMR  relaxation  rate,  which  averages  equally 
over  all  wave  vectors,  is  expected  to  disappear  or  be 
strongly  reduced.  It  was  shown  in  Ref.  [15]  that 
Josephson  currents  from  FSs  with  different  signs  of  the 
OP  interfere  destructively,  and  the  net  phase  corresponds  to 
the  sign  of  the  FS  with  the  higher  normal  conductance.  In 
the  constant  relaxation  time  approximation,  both  in-plane 
and  out-of-plane  conductivities  are  dominated  by  the  elec¬ 
tron  pockets.  This  is  unfortunate,  since  there  would  other¬ 
wise  be  a  77  phase  shift  between  the  ab  and  c  tunneling, 
and  corner-junction  experiments  could  be  used,  as  in 
high- 7/  cuprates. 

In  summary,  we  argue  that  the  fermiology  found  in 
doped  LaAsFeO  gives  rise  to  strong  but  broad  antiferro¬ 
magnetic  spin  fluctuations  near  the  M  point  in  the  Brillouin 
zone,  while  the  tendency  to  magnetism  existing  at  zero 
doping  is  suppressed.  These  fluctuations,  while  too  broad 
to  induce  a  magnetic  instability,  are  instrumental  in  creat¬ 
ing  a  superconducting  state  with  OPs  of  opposite  signs  on 
the  electron  and  hole  pockets. 
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